Spin- Transfer Torques in Helimagnets 
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We theoretically investigate current-induced magnetization dynamics in chiral-lattice helimagnets. 
Spin-orbit coupling in non-centrosymmetric crystals induces a reactive spin-transfer torque that has 
not been previously considered. We demonstrate how the torque is governed by the crystal symmetry 
and acts as an effective magnetic field along the current direction in cubic B20-type crystals. The 
effects of the new torque are computed for current-induced dynamics of spin-spirals and the Doppler 
shift of spin waves. In current-induced spin-spiral motion, the new torque tilts the spiral structure. 
The spin waves of the spiral structure are initially displaced by the new torque, while the dispersion 
relation is unaffected. 



I. INTRODUCTION 

Current-induced magnetization dynamics continue to 
be a very active research area due to potential appli- 
cations in future electronic devices. In metallic ferro- 
magnets, the magnetization can be manipulated via the 
spin-transfer torque (STT), which arises due to a mis- 
alignment between the spin polarization of the current 
and the local magnetization directional 2 - Slonczewski and 
Berger were the first to predict the existence of the STT 
effect which was later demonstrated in several exper- 
iments. 1 ' 2 The anticipated application potential of the 
STT effect lies in the development of electro-magnetic 
devices that utilize a current-induced torque instead of 
external magnetic fields to manipulate the magnetiza- 
tion. 2 

The magnetization dynamics of an itinerant fer- 
romagnet is described by the Landau-Lifshitz-Gilbert 
(LLG) equation extended to include the current-induced 
torques) 2 ^ 



rh = —7m x H e ff + am x m + t. 



(1) 



Here, m = M/M s (M s = |M|) is the unit direction vec- 
tor of the magnetization M, H e ff = —SF/SM. is the ef- 
fective field found by varying the free energy F\M] with 
respect to the magnetization, a is the Gilbert damping 
coefficient, 7 is (minus) the gyromagnetic ratio, and r 
describes the current-induced torques. In the absence of 
intrinsic spin-orbit coupling (SOC), the torque becomes 



= -(l-/3mx)(v s - V)m. 



(2) 



In Eq. ([2]), the first term is the adiabatic torque, while the 
second term (parameterized by (3) is the non-adiabatic 
torque. The vector v s is proportional to the current 
density, S, and its polarization, P: v s = —TiP J /2bsq. 
Here, sq is the total equilibrium spin density along — m, 
and e is the electron charge. The torque in Eq. ([2]) 
treats the ferromagnet within the exchange approxima- 
tion, which assumes that the exchange forces only depend 
on the relative orientation of the spins. This assumption 
is believed to be valid in metallic ferromagnets, including 
disordered systems in which impurities couple to the spin 



degree of freedom through random magnetic moments or 
spin-orbit coupling. In this case, impurity averaging re- 
stores the spin-rotational symmetry of the system. Re- 
cently, in systems with a broken spatial inversion sym- 
metry, the intrinsic SOC in combination with an external 
electric field have been observed to induce an additional 
torque, such that r = t cx + t so &~— In general, the SOC- 
induced torque, t so , can be written as&2J£ 



-7m x H s 



(3) 



where the SOC field, H so , is proportional to the electric 
field and its orientation is determined by the symmetry of 
the underlying crystal lattice and the direction of the ex- 
ternal electric field. Therefore, in contrast to t ox , which 
vanishes in a homogeneous ferromagnet, r so is finite even 
in this case. Several experiments have demonstrated that 
the SOC torque plays an important role in magnetization 
dynamics^ The underlying physics of the torque is that 
the SOC effectively acts as a magnetic field on the spins 
of the itinerant quasi-particles when an electric field is ap- 
plied to the system. The effective magnetic field induces 
an out-of-equilibrium spin density that yields a torque on 
the magnetization^^ 

In chiral magnets, the exchange interaction also con- 
tains an anisotropic term known as the Dzyaloshinskii- 
Moriya (DM) interactionj 11 ' 12 The DM interaction arises 
due to the characteristic crystalline asymmetry of the 
chiral magnet in combination with the SOC, and in cu- 
bic B20-type crystals, it leads to the formation of a spin 
spiral in the magnetic ground state. We refer to these 
systems as helimagnets. Helimagnets have recently at- 
tracted substantial interest because topological nontriv- 
ial spin structures, known as skyrmions, have been ob- 
served in such systems under the application of weak ex- 
ternal magnetic fieldsJ^r— Current-induced responses of 
the formed skyrmion lattice to current densities that are 
over five orders of magnitude smaller than those typically 
observed in conventional ferromagnetic metals have re- 
cently been observed experimentally— i2i To understand 
this striking feature of helimagnets, numerical simula- 
tions and a collective coordinate description have been 
applied to study the current-induced dynamics of spin- 
spirals and skyrmion lattices^ - — However, an impor- 
tant aspect of helimagnets is the absence of spatial inver- 
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sion symmetry, which implies that the magnetization ex- 
periences a SOC-induced torque given by Eq. None 
of the recent theoretical works have included or derived 
the explicit form of the term in Eq. (J3J. Its particular 
form and importance for the current-induced dynamics 
of hclimagnets therefore remain unknown. 

In the present paper, we derive the form of the torque 
in Eq. ([3]) for cubic non-centrosymmetric (B20-type) 
compounds. An important example of such a system 
is the chiral itinerant-electron magnet MnSi, which was 
the first system in which a two-dimensional lattice of 
skyrmions was observed. The effects of the SOC torque 
are studied for two different cases: current-induced spin- 
spiral dynamics and the Doppler shift of magnons that 
propagates along the spiral structure. We observe that 
for current-induced spin-spiral motion, the new torque 
yields an enhanced tilting of the spiral structure, while 
the torque does not affect the Doppler shift of spin waves 
except to induce an initial translation of the spiral struc- 
ture. We also briefly discuss the effect of the SOC torque 
on the skyrmion lattice dynamics. 

This paper is organized in the following manner. SecHIl 
provides a derivation of the SOC torque in Eq. for 
cubic B20-type crystals. Sec. |TTT] discusses the effects of 
the SOC torque on current-induced spin-spiral motion 
and the Doppler shift of spin waves that propagates along 
the spin-spiral. We conclude and summarize our results 
in Sec. El 



II. DERIVATION OF THE SOC TORQUE 

In deriving the explicit form of the torque in Eq. ([3]) , 
we are guided by the Onsager reciprocity relations and 
Neumann's principle. Consider a system described by the 
parameters {qi\i = 1, . . . , N} for which the rate of change 
qt is induced by the thermodynamic forces fa = —dF/dqi, 
where F(qi, ?jv) is the free energy. Onsager's the- 
orem states that the response coefficients in the equa- 
tions <ji = Y^j=iLijfj are related by Ljj(H,M) = 
ei6jLji(— H, — M), where ti — 1 (e^ = —1) if qi is even 
(odd) under time reversal^ m and H represent any pos- 
sible equilibrium magnetic order and an external mag- 
netic field, respectively. In the present paper, the re- 
sponses of the itinerant ferromagnet are described by the 
time derivative of the unit vector along the magnetiza- 
tion direction, rh, and the charge current density, J ' . 
The associated thermodynamic forces are the effective 
field scaled with the magnetization, f m = M s H c ff, and 
the electric field, f q = E, respectively, and the equations 
describing the dynamics in the linear response regime are 
determined by: 

(j) = ( L c L c) (t) • (4 » 

The Onsager reciprocity relations imply that L miqj (m) = 
—L qjmi (— m). In addition to the symmetry requirements 



imposed by the reciprocity relations, the symmetry of the 
underlying lattice structure also decreases the number of 
independent tensor components. This fact is expressed 
by Neumann's principle, which states that a tensor repre- 
senting any physical property should be invariant with re- 
spect to every symmetry operation of the crystal's point 
group £L 

According to Eq. ([J}, the effect reciprocal to the adi- 
abatic and non-adiabatic torque in Eq. ^ is a charge 
current density induced by a time-dependent magnetic 
texture. To the lowest order in the texture gradients 
and the precession frequency, the induced charge current 
density in the exchange approximation is 26 

j * = ^ p r x ^-^J- m - (5) 

Here, e is the electron charge, P is the spin polarization 
of the current, a is the conductivity, and r,; is component 
i of the spatial vector. Because the exchange approxi- 
mation neglects any coupling (via intrinsic SOC) of the 
spins to the crystal structure, the above expression is 
fully spin-rotational symmetric and a textured magneti- 
zation, i.e., dm/dri ^ 0, is required to have a coupling 
between the momentum of the itinerant quasiparticles 
and the magnetization. If the effects of intrinsic SOC are 
considered, additional terms are allowed by symmetry in 
the phenomenological expansion for the pumped current. 
In particular, for inversion symmetry-breaking SOC, a 
homogenous magnetization pumps a charge current. To 
the lowest order in SOC and precession frequency, the 
expression for the pumped current then becomes 

jr ap =T hi m j +jr. (6) 

The second-rank tensor r\ij is an axial tensor because the 
current is a polar vector while the magnetization is an ax- 
ial vector, rjij is linear in the SOC coupling and vanishes 
in systems with spatial inversion symmetry. According 
to Neumann's principle, the particular form of rjij is gov- 
erned by the crystal structure and is determined by the 
following set of equations produced by the generating 
matrices [Rij] of the crystal's point group^ 

Vij I ^l-^-inRjmVnm • (7) 

Here, \R\ is the determinant of the matrix [Rij]. 

Let us now consider a cubic B20-type crystal. Its crys- 
tal structure belongs to the non-centrosymmetric space 
group P2i3, which has the cubic point group T . Common 
examples of cubic B20-type chiral magnets are MnSi, 
FeGe, and (Fe,Co)Si. From the symmetry relations in 
Eq. ([7]) , one then finds that the tensor 77^ is proportional 
to the unit matrix^ 

Vij = ?7<%, (8) 

where is the Kronecker delta. The tensor is isotropic 
because the high symmetry of the cubic crystal reduces 
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the number of independent tensor coefficients to the sin- 
gle parameter rj. Substituting this tensor into Eq. © 
and expressing the time derivative of the magnetization 
in terms of the effective field by applying the first term on 
the right-hand side of Eq. (JTJ) , one obtains the response 
matrix! 2 ^ 



L 



<\, III , 



j ex 

qiirij " 



(9) 



Here L cx 



are the response coefficients describing the 
process reciprocal to the STT in Eq. @, which have been 
previously derived in Ref. [26l The term proportional to 
r\ describes the process reciprocal to the SOC-induced 
torque in Eq. ([3]). Using the Onsager reciprocity rela- 
tions, we find that the SOC field takes the following form: 



H s 



(10) 



where ^ so = (2r]eso)/(fiaPM s ). Thus, the torque in- 
duced by the SOC in non-centrosymmctric cubic mag- 
nets acts as an effective magnetic field along the current 
direction. Note that the torque is reactive because it 
does not break the time reversal symmetry of the LLG 
equation. 



III. RESULTS AND DISCUSSION 

In the following, we investigate the effects of the 
SOC torque on current-driven spin-spiral motion and the 
Doppler shift of spin waves. Additionally, a brief discus- 
sion of how we expect the torque to affect the skyrmion 
crystal dynamics is presented. 

A. Spin-Spiral Motion 

To the lowest order in the magnetic texture gradients, 
the free energy density of a ferromagnet with broken spa- 
tial inversion symmetry can be written phenomenologi- 
cally as)2£ 



. Jij dm dm 



D 



di 



dr k 



(11) 



Here, Jij is the spin stiffness describing the exchange 
interaction between neighboring magnetic moments, and 
the term proportional to Dyt is the DM interaction. In 
Eq. (jllj) (and in what follows), summation over repeated 
indices is implied. The explicit form of the tensors Jy 
and Dijk is determined by the crystal symmetry. 

In cubic B20-type ferromagnets, the free energy den- 
sity becomes 



_, . J 9m 9m 
^ (m) = 2 ^7 ' d^ 



Dm • (V x m) 



(12) 



The free energy of the system, F [m] = J drJ 7 , is mini- 
mized by a helical magnetic order, where the wave vector 



of the spiral structure is determined by the ratio between 
the DM parameter and the spin stiffness: k = D/J. For a 
k- vector that points along the z axis, the magnetic order 
of the ground state is 



mo(z) = cos(fcz)x + sin(fcz)y, 



(13) 



where x and y are the unit direction vectors along the x 
and y axis, respectively. 

The action functional, 5[m], and the dissipation func- 
tional, i?[m], of the system are written a a 29 i 30 



S[m] 



dtdr Ai (rhi + v s • Vra,) 



J_ 

M, 



T(m) 



R[m] = J 



dtdr 



7m • H 



111 



v, • Vm 



(14) 
(15) 



Here, A(m) is the Berry phase vector potential of a mag- 
netic monopole, which satisfies tijkdAu/ drrij = rrii (e^-fe 
is the Levi-Civita tensor). The LLG equation in Eq. (flj, 
with t = t cx + t so , is determined by 



(16) 



8S __ 5R 
Sm 5m 

A previous study on spin-spiral motion demonstrated 
that the response of the structure to an applied current 
(along z) can be described by the tilting angle, £, and 
drift velocity, £, of the spiral structure.— To find an ap- 
proximate solution of Eq. (|16l) , we therefore employ the 
following variational ansatz: 

m(z, t) = cos(£(i))m (z - C(t)) + sin(£(t))4. (17) 



Substitution of this ansatz into Eq. (IT4|) and (|15l) and 
integration over the spatial coordinates yield an effective 
action and dissipation functional for the variational pa- 
rameters £(i) and ((t): 



S[C,€] 



R[C,t] 



dt 

; (J 



1 , Oe 
= ldt 2 



(C~«.) ksm£+ (18) 
k 2 cos 2 ^ - Dk cos 2 £ ) - jH so sin^, 



-v s k - k( 
a 



(19) 



The equations of motion for the variational parameters 
are 



ss[c,£} 

5C 



SR[U] 6S[C,£] 5R[(,i} 



(20) 



st K ^ 

We are interested in the steady-state regime in which £ 
approaches a constant value. In this regime, the drift 
velocity and the tilting angle are 



c 

sin(0 





a 
M 



, , 1 «s #si 

7 Jk - 2D \ \ a I k 



(21) 



.(22) 
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The expression for the drift velocity £ agrees with the ex- 
pression derived in Ref. [22|. The SOC torque does not af- 
fect the drift velocity because the SOC torque effectively 
acts similarly to the adiabatic torque, as can be observed 
from the expression for the action S[C, £] in Eq. (fl8l) . The 
adiabatic and SOC torques initiate a spiral motion when 
a current is applied. However, the motion is damped due 
to the intrinsic pinning effect caused by Gilbert damping 
in combination with the DM interaction. Thus, similar 
to what is observed for domain walls in conventional fer- 
romagnets, a non-adiabatic torque is required to obtain 
a steady-state spiral motion. An observable effect of the 
SOC torque is the modification of the tilting angle ob- 
served in Ref. [H by an amount of —M s H so / (Jk 2 — 2Dk). 



applied along the z axis): 



B. Doppler Shift of Spin Waves 



In ferromagnets with a homogeneous magnetization, a 
Doppler shift in the spin wave dispersion relation under 
the application of a current has been observed^ The 
frequency, w, of the spin wave is shifted by v s • q, where q 
is the wave vector of the magnon: ui = (7 J/Af s )q 2 +v s -q. 

Theoretical works on Goldstone modes in helimag- 
nets with a spin spiral predict that these modes are 
much more complicated than those in homogeneous fer- 
romagnets^ We refer to these Goldstone modes as heli- 
magnons. The dispersion relation of the helimagnons is 
highly anisotropic, with a linear wave- vector dependency 
parallel to the spin-spiral direction and a quadratic de- 
pendency in the transverse direction (in the long wave 
length limit). That is, the soft modes behave like anti- 
ferromagnetic magnons along the spiral, while ferromag- 
netic behavior is observed for modes propagating in the 
transverse plane. Thus far, no works have studied the 
effect of an applied current on the dispersion relation of 
helimagnons. 

To derive an effective action for the Goldstone modes, 
we describe the local fluctuations by £ and ( in Eq. (fT?) 
by allowing the parameters to be both position- and time- 
dependent: £ = £(r,£) and £ — A similar param- 
eterization was performed in Refs. 13211331 in the analysis 
of helimagnons. The parameter ( describes a local twist 
(around the z axis) of the spiral structure, while £ de- 
scribes a local tilting along the z axis. In the analysis of 
the Doppler shift, we neglect dissipation and disregard 
the dissipation function. Ref. [33j demonstrated that sim- 
ple closed-form solutions for the variational parameters 
can only be obtained for modes propagating along the 
spin-spiral direction. For simplicity, we therefore restrict 
our study to Goldstone modes that propagate along the 
z axis. Expanding Eq. (fT4")l to second order in £(z,t) 
and £(2, t), we obtain the effective action (the current is 



dtdz k( ( C + v s — — — v. 

az 



2M, 



dz 



k 2 C 



(23) 



- 7#so£. 



The equations of motion are obtained by varying the ac- 
tion with respect to £ and £, i.e., 5S/5( — 6S/6£ = 0, 
which results in two coupled equations for the variational 
parameters: 



C(z,t)+v. 



i(z,t)+v t 



dgz,t) 

dz 



dz 



kM s 
, 7 



k 

~/Jk d 2 C(z,t) 
' M, dz 2 



(24) 
(25) 



Let us first consider the homogenous part of the equa- 
tions and neglect the two last terms on the right-hand 
side in Eq. (|24[) . Substitution of a plane wave ansatz of 
the form [Co £o] T exp(i(qz — uit)) into the equations yields 
the following dispersion relation: 



7J 
M, 



Vk 2 



v s q. 



(26) 



We see that the STT results in a Doppler shift similar 
to what is observed for spin waves in conventional ferro- 
magnets. In the long wave length limit, q — > 0, a linear 
dispersion relation is obtained: cj = (p/J /M s )kq + v s q. 
The SOC-induced torque only appears as a source term 
in the nonhomogeneous equations. The particular solu- 
tion (PS) of the nonhomogeneous equations in Eq. (p4|) 
and (I25l) is 



COM) 



.PS 



{l/k)H so )t 




(27) 



This solution describes a displacement of the spiral struc- 
ture induced by the adiabatic and SOC torques. How- 
ever, this current-driven spin-spiral motion is damped 
when dissipation is considered due to the intrinsic pin- 
ning effect. Thus, the SOC torque (together with the 
adiabatic torque) only causes an initial translation of the 
spin-spiral. 



C. Skyrmion Crystal Dynamics 

In helimagnetic thin-film systems, skyrmions have 
been observed under the application of a weak exter- 
nal magnetic field, B, perpendicular to the thin film. 
Each skyrmion has a vortex-like magnetic configuration, 
where the magnetic moment at the core of the vortex 
is anti-parallel to the applied field while the peripheral 
magnetic moments are parallel. From the peripheral mo- 
ments to the core, the magnetic moments swirl up in 
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a counter-clockwise or clockwise manner. The formed 
skyrmions arrange themselves in a crystalline structure, 
a two-dimensional skyrmion crystal. 

Recent experiments have revealed curren-driven 
skyrmion crystal motion at ultra low current densities. 20 
The motion of a skyrmion lattice is only weakly affected 
by pinning, which is in stark contrast to observations for 
current-induced domain wall dynamics in conventional 
ferromagnets. A theoretical work has indicated that the 
pinning-free motion arises because the skyrmion lattice 
rotates and deforms to avoid the impurities. 2 - However, 
all analyses of current-driven skyrmion crystal motion 
have disregarded the SOC torque. 

Sec. ITlI showed that the SOC torque acts as an effective 
field along the current direction. For a current applied 
along any direction in the thin film, the expected conse- 
quence of the SOC torque is therefore that this torque 
leads to a small correction to the external magnetic field 
that stabilizes the two-dimensional skyrmion lattice such 
that the total field becomes Ht = B+H so . The expected 
response of the skyrmion crystal to this perturbation is 
a rotation of the 2D lattice structure that aligns the core 
magnetic moments anti-parallel to Ht- To confirm our 
predictions, a more thorough numerical simulation of the 
magnetic system is required, which is beyond the scope 
of the present work. 

IV. SUMMARY 

In this paper, we performed a theoretical study of 
current-induced torques in cubic non-centrosymmetric 



helimagnets. We demonstrated that due to the broken 
spatial inversion symmetry, the SOC induces a reactive 
magnetization torque that has not been considered in 
previous studies of current-induced dynamics in helimag- 
nets. The specific form of the SOC torque is determined 
by the symmetry of the underlying crystal lattice and 
acts as an effective magnetic field along the current di- 
rection in B20-type chiral magnets. 

The consequences of the SOC torque are studied for 
two different cases: current-induced spin-spiral motion 
and the Doppler shift of helimagnons. During the 
current-driven spin-spiral motion, the SOC torque yields 
an enhanced tilting of the spin-spiral structure, while the 
velocity is not affected. The dispersion relation of a heli- 
magnon that propagates along the axis of the spin-spiral 
is not affected by the SOC torque except to induce an 
initial translation of the spiral structure. 
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